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Abstract 

Hydrodynamic flows in biological systems are often generated by active chiral processes near or 
on surfaces. Important examples are beating cilia, force generation in actomyosin networks, and 
motile bacteria interacting with surfaces. Here we develop a coarse grained description of active 
chiral films that captures generic features of flow and rotation patterns driven by chiral motors. 
We discuss force and torque balances within the film and on the surface and highlight the role of 
the intrinsic rotation field. We arrive at a two dimensional effective theory and discuss our results 
in the context of ciliary carpets and thin films of bacterial suspensions. 
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The building blocks of biological systems are chiral molecules such as proteins and DNA. 
The emergent collective behaviors of these molecular components give rise to active dynamic 
processes in cells and tissues which can also reflect these molecular chiralities. A key example 
is the breaking^ of left-right symmetry in organisms with a well-deflned handedness during 
development jlMsl] (i.e. the heart is on the left side in humans). Many pattern forming 
events in biological systems are governed by active chiral processes. For example, it has 
been shown that the rotating beat of cilia, which drives chiral hydrodynamic flows, is at 
the basis of left-right symmetry breaking in vertebrate animals (#-1611. In addition, active 
chiral processes have been observed in the cytoskeleton of cells [7|, [sl]. The cytoskeleton is 
a gel-like network of elastic chiral fllaments and other components such as motor proteins. 

n n 

The interactions of motors and fllaments drive movements and intracellular flows [9|, I10| . 
which can have chiral asymmetries j?. 



11 



Active chiral processes are typically observed on surfaces or at interfaces. Key examples 
are carpets of beatin g ci lia driving hydrodynamic flows parallel to the surface, on which 
they are attached U^, ll5fl, and rotating motors on a surface 



Many different beating 



patterns of cilia exist, which in general are chiral and often exhibit rotating movements, as 



is the case for the cilia that govern the left-right symmetry breaking of organisms p, ll9|| . 



Many microorganisms possess carpets of cilia on their outer surface, which are used for 



self-propulsion along helical trajectories in a fluid p^. Recently, the collective behavior of 
swimming E. coli bacteria on solid surfaces was studied These bacteria possess rotary 
motors, which drive the rotation of helical flagella relative to the cell body. This relative 
rotation provides a motor for self-propulsion in a fluid. Due to conservation of angular 
momentum, cell body and flagella exert equal but opposite torques on the surrounding 
fluid, see Fig. [I](a). Thus, the bacteria act as torque dipoles. Such flagellated E. coli were 
placed on a solid agar surface j2o|, where they produced a thin fluid layer that promotes 
their motility. Some bacteria attach to the surface while others swim in the fllm. Those 
bacteria which were attached to the surface also exerted torques on the solid substrate. 
Large scale chiral flow patterns were reported as a result of bacterial activity The 
bacterial suspension therefore represents an internally driven active chiral fllm supported by 
a surface. 

Active fluids and gels have been studied in the framework of hydrodynamic theories both 
as a paradigm for the cell cytoskeleton and for suspensions of active swimmers |2l|-|26|. 



2 



Such approaches are based on hquid crystal hydrodynamics (27l-l32| driven out equihbrium 
by internal active processes. It has been shown that stresses generated by active processes 



can give rise to a rich variety of dynamic patterns and flows 



38. 
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10. 133l-l37|. Similar approaches 



Chiral effects have so far 



40| . A study of active chiral 



have also been used for the study of granular systems 
been discussed mainly in passive systems and in granular gases 
processes in fluids and gels is so far lacking. 

Here, we develop a generic theory for active chiral fllms. We consider a thin fllm of a 
suspension of chiral motors on a solid surface. Our work is based on a general discussion 



of the bulk properties of active chiral fluids 



4l|. Active chiral processes are introduced as 



torque dipoles arising from counterrotating objects, such as the cell body and the flagella of 
bacteria, see Fig. [l](a). 

We start by discussing the bulk properties of a fluid with mass density p and center 
of mass velocity v in which active processes take place. Linear momentum and angular 
momentum conservation can be expressed as 



(1) 
(2) 



where Einstein's summation convention is implied. Here pv^ is the momentum density. The 
density of angular momentum /^^j is described by an antisymmetric second rank tensor. 
Externally applied force and torque densities are denoted by and r^^^, respectively, is 
a position vector, is the total stress and M^^^^ is the total angular momentum flux. The 
total stress a^^J = —PSa/s + 0-^/3 + describes momentum fluxes and can be decomposed 
in the isotropic pressure P, the symmetric traceless stress and the antisymmetric stress 



a^^. The total angular momentum density consists of an orbital part 



r^Va) and a 



spin angular momentum density l^^ = /^^j — p{raVp — T^Va)^ see |27|, l28|, [321, l40|. Similarly, 
the total angular momentum flux M^^^ = M^isj + {^acrf^ — "i^/s^^aj) ^he sum of fluxes M^/Sj 
and r^crf^ — "^/s^^aj ^P^^ orbital angular momentum, respectively. Note that the spin 
angular momentum density la/s and the spin angular momentum flux M^isj do not depend 
on the choice of coordinate system. We deflne an effective rate of intrinsic rotation Qa(3 of 
local volume elements via la(3 = la/B^s^jS^ where Ia(3j6 is the moment of inertia tensor per 



unit volume 



. The balance of the local torque then reads 



ext 
a(3 



(3) 




FIG. 1: (a) Schematic representation of a swimming bacterium with rotating flageha. (b) Chiral 
motor consisting of two counterrotating spheres at distance d. Arrows indicate senses of rotation. 

Eq. (|3j) shows that spin angular momentum is not conserved but can be exchanged with 
orbital angular momentum by antisymmetric stress. 

Active chiral processes introduce contributions to the total stress and to fluxes of spin 
angular momentum. We discuss the effect of an active chiral process by considering a torque 
dipole Tai3{r) in a fluid, located at r = 0, see Fig. [I](b). This torque dipole is built from two 
torque monopoles iqCa/s^Pj^ of strength q separated by a small distance d in the direction 
of the torque axis given by the unit vector p: 



Note that ra(3 is invariant under the transformation p — p, which implies a nematic 
character. 

Active internal torques and forces are introduced similarly to external ones. Therefore the 
torque dipole r^/^ can be interpreted as an active contribution M^^^^ = —qdea/ssPsPj^i'^) to the 
spin angular momentum flux, which obeys d^M^^^^ = Ta(3- In a suspension of many identical 
torque dipoles at positions r*^^^ and with orientations p^^^ the active angular momentum 



- -qdeaf3uPuPjd^S{r) . 



(4) 



fluxes are 




act 
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- C^a(3SP6Pj + C^a^j , (5) 

where ( = —Sqdn{r) and = {S — l)qdn{r)/3 describe the strength of a nematic and an 
isotropic active contribution to Ma(3j^ respectively, and n(r) is the dipole density. Here and 
below p, with = 1, denotes a coarse grained nematic director and is a nematic order 
parameter which obeys S{paPj± — {l/3)6a(3) =< P^aP^^ > —{l/3)Sa(3^ where the average is 



taken over a volume element 



30]. 



The constitutive equation for the spin angular momentum flux then has the form 

Ma(3j = l^d^^a(3 + C^a/BSPsPj + C^a(3j • (6) 

Here the term proportional to the phenomenological coefficient k describes passive dissipative 
processes. Furthermore we ignore other passive couplings {32, 41|. Here and below we have 
neglected inert ial terms. We complement Eq. (|6]) by the constitutive equations for the 
stresses of a passive incompressible fluid: 

<^ = 2r]\Q^^ - uj^^) , (8) 

where Ua/s is the traceless part of the strain rate Ua/s = {d^vi^ + d^v^) /2 and 00^(3 = 
{daVi3 — df^Va) /2. Here, rj is the shear viscosity and 77' is a rotational viscosity. The pressure 
P plays the role of a Lagrange multiplier and is imposed by the incompressibility condition 
djVj = 0. Other coupling terms, known to exist in liquid crystals, have been neglected for 
simplicity. 

We now study the stresses and flows generated by active chiral processes in a thin fluid 
fllm of height /i on a solid substrate, see Fig. [2l Note that the fllm is not symmetric with 
respect to z ^ h — z because the two surfaces at z = and z = h differ. Integrating the 
force balance Eq. ([T]) we obtain 

h 

0=11 dzdpalf = d,a^ + y:ro , (9) 



where the bar denotes an average over the fflm height, a^^^ = (l/h) dz(j\^^ ^ and = 
(jlf{z = h) — (jlf{z = 0), where i = x^y and j = see Fig. [2l The shear stress at the 
film surfaces is written as 

l'jiX = -^tV^ + Ca^i. , (10) 
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which describes friction forces due to relative motion with the substrate by the coefficient 
and due to relative rotation with the substrate by the coefficient ^q. From the torque 
balance Eq. ([3]) and ignoring inertial terms, we find 

= -2a^^^ + d,M^0, + ^M^0X ■ (11) 
The torque densities on the surfaces are 

-^Ma(3z\o = —^r^a(3 + Cs^a(36pzp6 + Cs^a(3z , (12) 

where frictional torques are described by the coefficient ^r- The torque exerted by active 
chiral processes on the boundaries is described by the coefficients Cs and C^. Averaging 
over the height of the film also leads to new coefficients for the active in-film angular mo- 
mentum ffuxes M(y,(3j = Cf^a(36p6pj + Cf^a(3j- Finally, the externally applied shear stress on 
the boundaries is a^^^ ^ = [r] -\- Ti')vi\^/h -\- 2r]^Cliz. In the absence of external shear 
stresses Vil^ = —2hr]'Qiz/{r] + ij'). Similarly, = P^^^ is the externally applied pressure 
and therefore P = —2r]diVi. We obtain equations for the fiow and rotation field, 

= dj [{3r] - r]')diVj + [r] + r]')djVi] + 277'^^% - it^i + in^iz , (13) 

= -{At]' + ir)^xy + 27y'(5^^2/ " ^2/^^) + ^^j + djQe^yzPzPj + Cs^xyzpl + Cs^xyz (14) 

^ = ~ ( + ) ^iz + l^d'^j^iz + djCfCizSpSpj + djCf^izj + Cs^izSpSpz , (15) 

where we used the thin film approximation diV^ <C d^Vi. 

Motivated by rotating cilia or bacteria attached to a surface, we now consider the case 
of a thin active film containing chiral motors that are on average aligned along the vector 
p = (cos sin 0, sin sin 0, cos ^). Here we have introduced the angle 9 which describes the 
average tilt with respect to the surface normal vector and the angle (j) which specifies the 
tilt direction with respect to the x-axis. We first consider a homogeneous distribution of 
motors in an infinite system where all spatial derivatives vanish. In this case Eq. (fT3l) reads 
Vi = {^n/^t)^iz- Using Eq. ^ we find 

^izjPjPz /-i^x 
' ~ ^r + ^VV'/iV + V') ' 

A fiow with a velocity |v| oc sin(2^) is generated. The fiow direction is perpendicular to the 
direction of tilt with respect to the surface normal vector. This case describes the collective 
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FIG. 2: Schematics of a thin fluid fllm of height h that contains chiral motors. Motors are torque 
dipoles that consist of counterrotating spheres (see Fig. [11(b)), one of which is attached to the 
surface. The other rotates as indicated by the arrows. The average motor direction is described 
by the vector p which is tilted at an angle 6 in ^-direction. 

generation of flow by carpets of cilia on a surface j^l, Q • Note that and therefore also 
the generated flow vanish in fllms that are symmetric with respect to z ^ h — z. Ciliary 
carpets do have this asymmetry and can therefore generate flows [6]. Furthermore, we flnd 
from Eq. f fT4l ) an intrinsic rotation rate Qxy = (CsCos^^ + C^)/(477' + ^r) generated by the 
rotating motors with vorticity CUxy = in the hydrodynamic flow fleld. 

As a second example we consider a circular patch of radius R that contains active motors. 
This is described by position dependent coeflScients Cs(i^) = Cs^{R — \^\) and by similar 
expressions for (^^, (f and where B(r) is the Heaviside function. We numerically solve 
Eqs. ( IT3l) . ( IT4l) and f fTSl) with periodic boundary conditions for boxsize L = 4i? using Fourier 
transforms, see Fig. [3l The velocity fleld for a tilt angle ^ = is displayed in Fig. [3](a). The 
patch generates a chiral flow fleld driven by the intrinsic rotations Cl^y. The stream lines 
are concentric circles. The flow velocity is maximal at the edge of the patch and decays 
exponentially outside. There is no net transport across the patch. If the motors are tilted 
with 9^0 along the y axis, a net transport in x direction across the patch with velocity 
proportional to sin(26^) appears in addition to the circular flow, see Fig. [3](b). 

We have developed a coarse grained description of flow patterns generated in thin active 
fllms in which chiral processes take place. We have shown that if the chiral motors are 
tilted with respect to the surface normal vector, directed flows can be generated over large 
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a) 



b) 




FIG. 3: Flow fields (vectors) generated by a circular patch of chiral motors (blue), a) Average 
motor axis p is perpendicular to the surface, i.e. 6 = 0. b) Average motor axis p tilted at 
angle = 7r/4 in ^/-direction. Flow fields were determined numerically in a square box of size L 
with periodic boundary conditions. Parameter values are rj' = r] = h^Hi = h'^^^t — —h^,Q = (^r? 
( = -hCs = IOC' = -10/iCs, L = 20h. 

distances in a direction perpendicular to the tilt direction. The velocity of this net flow 
is maximal for tilt angles of 6^ = 45°. This result accounts for the generation of net flows 



and left-right symmetry breakin, 



Kupffer's vesicle in zebrafish 



g by carpets of cilia in the mouse ventral node 



42 



and 



]. Note that in both cases a tilt angle with respect to the 
surfa^ normal can be defined il9|]. Interestingly, reported tilt angles vary between 30° and 
50° [43], which is close to the angle of maximum transport velocity. Moreover in the case 
of finite patches of chiral motors, intrinsic rotation rates drive chiral flows along the edge of 
the patch. This phenomenon was reported in recent experiments on bacterial films on solid 
surfaces j2o| . 

Our theory highlights the role of the intrinsic rotation field Q^/s for active chiral pro- 
cesses. In particular we find net flows without vorticity generated by intrinsic rotations in 
a homogeneous system, see Eq. f fT6l) . Note that in passive bulk fluids Qa/3 converges to the 
flow vorticity after a relaxation time that in general is short [32]. In fllms of chiral motors, 
however, the intrinsic rotation rate differs from the vorticity even at steady state, and can 
create effects near interfaces and surfaces. Our work shows that the hydrodynamic flow 
velocity v and the intrinsic rotation rate Qa(3 in active chiral fllms are coupled by the ac- 
tive processes and by the boundary conditions. This can give rise to complex flow patterns 
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generated by carpets of chiral motors. 

Our theory permits a simplified description of the flows generated by ciliary carpets. It 
will be interesting to expand this theory to include phase relationships between adjacent 
motors to capture the self-organization of beating cilia in metachronal waves fisl . 
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